plot[(((((ceiling[(x!)^0.5]+k)^2-x!)^0.5-ceiling[((ceiling[(x!)^0.5]+k)^2-x!)^0.5])^(1/((10)^100+1)))^2)] for k=0